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the free parameters of this model are obtained. 
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1 Introduction 

The attempt to unify all fundamental interactions resulted in the devel- 
opment of multidimensional theories like string derived field theories ^ El 
El EJ 03 EJ, related brane- world theories El El El) and (related or not) 
Kaluza-Klein theories [111 1121 [T31 1141 115j . Among these theories, there are 
some in which the gauge coupling constants may vary over cosmological time 
scales. On the other hand, a theoretical framework based on first principles, 
was developed by Bekenstein |16j and later improved by Barrow, Sandvik 
and Magueijo in order to study the possible time variation of the fine 
structure constant. Furthermore, this model was generalized to study the 
variation of the strong coupling constant by Chamoun et al.|18j. 

Different versions of the theories mentioned above predict different time 
behaviours of the gauge coupling constants. Thus, bounds obtained from 
astronomical and geophysical data are an important tool to test the validity 
of these theories. 

The experimental research can be grouped into astronomical and local 
methods. The latter ones include geophysical methods such as the natural 
nuclear reactor that operated about 1.8 10 9 years ago in Oklo, Gabon ^SlEOl 
121) , the analysis of natural long-lived j3 decayers in geological minerals and 
meteorites |22| I2H ( 12*1] and laboratory measurements such as comparisons of 
rates between clocks with different atomic number UJ3 HH1 E7J EE1 EH1 E0]. 
The astronomical methods are based mainly in the analysis of spectra form 
high-redshift quasar absorption systems j3Tl 1321 1331 151 1551 I5n1 1571 I5H1 1551 HOI 
141) . Although, most of the previous mentioned experimental data gave null 
results, evidence of time variation of the fine structure constant was reported 
recently from high-redshift quasar absorption systems [3J3 EH E3 EfH EH 
I4U| . However, other recent independent analysis of similar data |411 1421 
1431 144j found no variation. On the other hand, measurements of molecular 
hydrogen |38| 140] reported a variation of the proton to electron mass /i = 
— . Furthermore, the time variation of the gauge coupling constants in the 
early universe can be constrained using data from the Cosmic Microwave 
Background (CMB) GJJ EH1 and the primordial abundances of light 

elements 031501121]. 

Damour and Polyakov [1] proposed a string dilaton model in which the 
variation of the gauge coupling constants is driven by the time evolution of 
the dilaton field. In this model, the dilaton remains massles and in conse- 
quence all coupling constants and masses of elementary particles are time 
dependent. In this chapter we limit ourselves to the variation of the fine 
structure constant. We compare the prediction of the dilaton model with 
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all available data on time variation of the fine structure constant but the 
nucleosynthesis data. In section [21 we describe carefully the data set consid- 
ered. In section |31 we briefly review the Damour and Polyakov proposal |H 
and obtain the analytic expression for the prediction of the time variation of 
the fine structure constant including the renormalization group corrections. 
Finally, in section |1] we present the results of the statistical analysis and 
briefly discuss our conclusions. 



2 Bounds from astronomical and geophysical data 

In this section, we review all availabe bounds on time variation of the fine 
structure constant. We discuss the relation between the observable quanti- 
ties and the variation of a in each case. We also describe very carefully the 
error we consider for each data. 



2.1 The Oklo Phenomenon 

One of the most stringent limits on time variation of the fine structure 
constant a follows from an analysis of isotope ratios in the natural uranium 
fission reactor that operated 1.8 x 10 9 yrs ago at the present day site of the 
Oklo mine in Gabon, Africa. From an analysis of nuclear and geochemical 
data, the operating conditions of the reactor could be reconstructed and the 
thermal neutron capture cross sections of several nuclear species measured. 
In particular, a shift in the lowest lying resonance level in 149 Sm : A = 
£,i49(Okio) _ ^,i49(now) ^ e derived from a shift in the neutron capture 
cross section of the same nucleus |20( I19| . The first estimate of a change in 
the resonance energy was peformed by Damour and Dyson ^H]- This bound 
was re-examined by Fujii et al [2U] using new samples of 149 Sm, 155 Gd and 
157 Gd. In their analysis they take the effect of contamination into account, 
assuming the same contamination parameter for all samples. They obtain 
the following bound: 

A = E? Ho - El ^ = (9 ± 11) 10' 3 eV (1) 

The shift in A can be translated ^HI into a bound on a possible difference 
between the value of a during the Oklo phenomenon and its value now, as 
follows: 

A = a^^ (2) 
da a 
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where = 10 . The dependence of the resonance energy with other 
fundamental constants has been analized by Damour and Dyson ^S] and 
Sisterna and Vucetich 52 . 



2.2 Long-lived (3 decayers 

The half-life of long-lived j3 decayers has been determined either in labora- 
tory measurements or by comparison with the age of meteorites, as found 
from a decay radioactivity analysis. In table 12.21 we show for three 
diferent decayers: 187 Re, 40 K, 87 Rb. Sisterna and Vucetich [2S] have derived 
a relation between the shift in the half-life of long lived (3 decayers and a 
possible variation between the values of the fundamental constants a, Aqcd 
and Gf at the age of the meteorites and their value now. In this chapter, 
we only consider a variation and therefore, the following equation holds: 

AA Aa 

— = a— 3 

A a 

where a = 21600, 46, 1070 for 187 Re, 40 K, 87 Rb respectively 



Table 1: The table shows the (3 decayer, the difference of half life measured in 
the laboratory and from meteorites, the corresponding error and reference. 



(3 decayer 


AX 

A 




Reference 


l87 i2e 


-1.6 10" 2 


1.6 10~ 2 


El 


m K 





1.3 10- 2 




87 RB 





1.3 10~ 2 


22. 



2.3 Atomic clocks 

The comparison of different atomic transition frequencies over time can be 
used to determine the present value of the temporal derivative of a. In- 
deed, the more stringent limits on the variation of a are obtained using this 
method. 

Hyperfine transition frequencies have the following dependence with a: 
VH yv ~ a 2 ——RooCFre^uZ) (4) 

where \x is the nuclear magnetic moment, /x^ is Bohr magneton, is Ry- 
dberg's constant, m p and m e are the proton and electron mass and Free 
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is the relativistic contribution to the energy. In such way, the comparison 
of rates between clocks based on hyperfine transitions in alkali atoms with 
different atomic number Z can be used to set bounds on a k where k 



depends on the frequencies measured and fj,A, t refers to the nuclear mag- 
netic moment of each atom. As explained above, we are considering only a 
variation, but it is important to keep in mind that this type of experiments 
actually constrain a combination of a and other fundamental quantities. 
The first three entries of table 12.31 show the bounds on — obtained com- 

1 1 a 

paring hyperfine transition frequencies in alkali atoms. On the other hand, 
an optical transition frequency has a different dependence on a: 



where B is a numerical constant assumed not to vary in time and Fi{a) 
is a dimensionless function of a that takes into account level shifts due to 
relativistic effects. Thus, comparing an optical transition frequency with 
an hiperfine transition frequency can be used to set bound on d k ^j^- 
Again, we will only consider a variation. Different authors |28l 1291 I3L)| have 
measured different optical transitions and set bounds on the variation of a 
using different methods. Fischer et al have considered the joint variation 
of a and We have reanalized the data of ref. [2JJ|, considering only a 
variation, yieding the fifth entrie of table 12.31 On the other hand, Peik et 
al, have measured an optical transition frequency in 171 Y b+ with a cesium 
atomic clock. Furthermore, they perform a linear regression analysis using 
this result toghether with other optical transition frequency measurements 
|28| I29j . On one hand, a linear regression analysis with three points has 
no statistical significance. On the other hand, we have already considered 
the other data. Therefore, we have also reanalized the data, using only the 
comparison between Yb+ and Cs frequency, yielding the sixth entrie in table 



2.4 Quasar absorption systems 

Quasar absorption systems present ideal laboratories to search for any tem- 
poral variation in the fundamental constants. The continuum spectrum of 
a quasar was formed at an epoch corresponding to the redshift z of main 
emission details specified by the relationship: 



MA 2 




(5) 



ESI 



Kbs — ^lab (1 + Z) . 



(6) 



Quasar spectra of high redshift show the absorption resonance lines of 
alkaline ions like CIV, Mgll, Fell, SilV and others, corresponding to the 
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Table 2: The table shows the clocks compared, the value of — and its 
corresponding error, the time interval for which the variation was measured 
and the reference. 



Frequencies 


Aa 

Q 






At(yr) 


Reference 


Hg+ and H maser 







1.4 10" 14 


0.38 


US] 


Cs and Rb 


8.4HT 


15 


13.8 10~ 15 


2 


M 


Cs and Rb 


-210" 


16 


810~ 15 


5 


nu 


Hg and Cs 







2.4 10~ 15 


2 


ES 


H and Cs 


5.710" 


15 


11.2 10~ 15 


5 


EH] 


Yb and Cs 


-1.610 


-15 


5.9 10" 15 


2.8 


M 



5*1/2 — ¥ P3/2 (-^l) an d 5*1/2 — * P1/2 (^2) transitions. The relative magnitude 
of the fine splitting of the corresponding resonance lines is proportional to 
the square of the fine structure constant a to lowest order in a. 

AA Ai - A 2 2 /r,\ 
— = ^—~a (7) 

Therefore, any change in the value of a at redshift z with respecto to 
the laboratory value, can be measured from the separation of the doublets 
A A as follows: 



Aa _ x 

a A 



( A )z _y 



1 AA \ 
v A / no 

where c r is a correction term which depends on the ion considered. This 
method is known in the literature as the Alkali Doublet (AD) method. Sev- 
eral authors have applied this method to SilV doublet absorption lines sys- 
tems at different redshifts (1 < z < 3.6). We show the average values they 
obtain in table |3] On the other hand, in order to perform our statistical 
analysis, we consider the individual data and corresponding errors for each 
absorption cloud which can be found in the references. 

Bahcall et al use strong nebular emission lines of O III to constrain 
the variation of a. Again, in this case, measuring the relative separation of 
the doublet give a constraint on a 2 at a given redshift. Assuming a linear 
variation of a with time, the authors find ^ = (—0.7 ± 1.4) x lO -4 . In 
this work , we want to test a theoretical model and thus the data must be 
model independent. Thus, we will consider the individual measurements on 
the variation of a belonging to the standard sample consisting in 42 quasar 
absorption systems, which are listed in table 4 of reference 141 j . 



Testing a string dilaton model with experimental and observational data 7 



Table 3: The table shows the redshift interval, the average value and stan- 
dard deviation of — in units of 10 _4 b obtained using the AD method and 
the corresponding reference. 



Zabs 


Aa 

a 

xlO" 4 


xlO" 4 


Reference 


2.6 < z < 3.6 


-0.35 


3.5 




2.8 < z < 3.05 


0.21 


1.4 


EH 


2 < z < 3 


-0.05 


0.13 


mm 


1.18 < z < 1.83 


-0.31 


0.85 





An improvement to the AD Method was proposed by Webb et al |331 153| . 
The new method called in the literature the Many Multiplet (MM) method 
compares transitions of different species, with widely differing atomic masses 
toghether with different transitions of the same species. If we consider a 
many electron atom or ion, the relativistic correction to the energy of the 
external electron can be written as: 



A~ (Z n a) 2 \E 



3/2 



1 



C(j,l) 



(8) 



J + 1/2 

where Z n is the nuclear charge, E is the electron's energy, j and I are 
the electrons's total and orbital angular momentum respectively. Moreover, 
C (j, I) is the contribution added due to the many body effects. Again, as 
in the case of the atomic clocks, the relativistic contribution is proportional 
de Zo? and therefore comparing different transitions of the same atom or 
transitions of different atoms, is a useful tool to put bounds on a. Moreover, 
the energy equation for a transition to the ground state, within a particular 
multiplet, at a redshift z reads: 



E K 



E c + 



+ K x L.S 



Z 



Oz\' 

a J 



1 



K 2 ( L.S) Zi 



-1 

a 



(9) 

where L and S are the electron total orbital angular and spin momentum 
respectively, E c is the energy of the configuration centre, Qi, K% y K2 are 
relativistic coefficients which have been accurately computed by Dzuba et 
al. |531 1541 155j . The limits over the fine structure constant variation are 
obtained fitting Voigt profiles to the absorption features in several different 
transitions. To the three usual fit parameters : column density, Doppler 
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width and redshift, they add In such way, this method makes possible 
to gain an order of order of magnitude in sensibility with respect to the 
AD method. On the other hand, the simultaneous fit of the Voigt profiles, 
can allow that unknown systematic errors hide under the variation of a. As 
mentioned before, this method provides the only results consistent with a 
time varying fine structure constant [331 EU ESI EH] • In a recent work [33] 
they estimate ^ = (-0.543 ± 0.116) x 10~ 5 for 128 absorption systems 
over the redshift range 0.2 < z < 3.7 confirming their previous results. For 
our analysis we consider the individual values obtained for each absorption 
system, listed in table 3 of reference EH[. As suggested by the authors we 
add to the each individual error listed in table 3 of ref. [SHI) an additional 
random error of 2.09 x 10~ 5 . 

On the other hand, the same method was applied by other authors |441 
I56| . who also use a stringent selection criteria of the samples, discarding 
weak and blended lines among others. They have obtained no variation of a 
for a high quality quasar spectra obtained using VLT over the redshift range 
0.4 < z < 2.3. Again, in this case, we use the individual estimates obtained 
for each redshift listed in table 3 of reference [SHj • On the other hand, the 
standard MM tecnique can be revised to avoid the deficiencies pointed out 
earlier and in the literature [JJJE7JIIS]- In fact, for ^ << 1 equation EH 
can be re-written as follows: 

Zi = Z a + K a Qi (10) 

where Z\ denotes the observed redshift and Qi the sensitivity coefficient 
corresponding to the lines i, and the slope parameter n a is given by: 

K a = -2{l + z a ) 11 

a 

In such way, it is possible to estimate the variation of a from linear 
regression analysis of the position of the line centroids in an absorption 
component. The accuracy of the regression analysis will be improved, if 
several absorption line samples are combined. This improved method called 
in the literature as Revised Many Multiplet (RMM) method, was applied by 
Levshakov [ST] and Quast et al [33] to a homogeneous sample of Fell lines 
at redshift z = 1.149 and z = 1.15 to obtain = (1-1 ± 1-1) x 10" 5 and 
4^ = (-0.4 ± 4.6) x 10~ 5 respectively. 

On the other hand, OH lines can provide precise constraints on cosmic 
evolution of a, the proton g factor and the ratio of electron to proton mass. 
Again, in this chapter, we will only consider a variation. Darling [SH] re- 
ported the detection of of the satellite 18 cm OH conjugate lines at 1612 
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and 1720 MHz in the z = 0.2467 molecular absorption system toward the 
radio source PKS 1413 + 135. Conjugate lines profiles guarantee that both 
lines originate in the same molecular gas. On the other hand, the 18 cm 
OH lines can be decomposed into a A-doubled term chich depends weakly 
on a and a hiperfine term which has a strong a 4 dependence. From these, 
sums and differences of lines can form pure A-doubled and pure hyper fine 
quantities: 



= v 1720 + 1/1612 = 2Aa - 4 
Au = 1/1720 - ^1612 = 2 (A + + A - ) a 4 

where A+ = 9.720 MHz and A" = 9.375 MHz and A = 11926MHz . In such 
way, Darling obtains ^ = (0.5 ± 1.3) x 10" 5 . 

Moreover, the ratio of the hyperfine 21 cm absorption transition of neu- 
tral hydrogen u a to an optical resonance transition is proportional to 
x = ® 2 9p^ where g p is the proton g factor. Thus, a change of this quantity 
will result in a difference in the redshift measured from 21 cm and optical 
absorption lines as follows: 

Ax _ Z opt - 2 2 1 



x (l + z) 



(12) 



So, combining the measurements of optical and radio redshift, a bound 
on a 2 can be obtained. Table 0] shows the bounds obtained by different 
authors using this method. This method has the inconvenience that it is 
difficult to determine if both radio and optical lines originate at the same 
absorption system. Thus, a difference in the velocity of the absorption clouds 
could hide in a variation of a. 

The ratio of the rotational transition frequencies of diatomic molecules 
such as CO to the 21 cm hyperfine transition in hydrogen is proportional to 
y = g p a 2 . Thus, any variation in y would be observed as a difference in the 
redshifts measured from 21 cm and molecular transition lines: 

Ay z mo i - z 2 i 



y (1 + 2) 



Murphy et al. |61j have placed upper limits at redshift z = 0.25, 
(-0.1 ± 0.22) x 10~ 5 and at redshift z = 0.68, ^ = (-0.08 ± 0.27) x 10" 5 . 
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Table 4: The table shows the absorption redshift, the value and standard 
deviation of ^ obtained comparing optical and radio lines in units of 10~ 4 
and the corresponding references 



Zabs 


Aa 
n 


V a ) 


Reference 


1.77 


-0.035 


0.055 


31. 


0.52 





0.6 


Em 


0.69 





1.4 


60. 



2.5 Cosmic Microwave Background 

Any variation of the fine structure constant a alters the physical conditions 
at recombination and therefore changes the cosmic microwave background 
(CMB) fluctuation spectrum. The dominant effect is a change in the redshift 
of recombination, due to a shift in the energy levels and, in particular the 
binding energy of Hidrogen. The Thompson scattering cross section is also 
changed for all particles, being proportinal to a 2 . A different value of a 
at recombination affects the CMB fluctuation spectrum in two ways: i) a 
shift in the Doppler peaks position and ii) a change in the amplitude of 
the Doppler peaks. On the other hand, the CMB fluctuation spectrum is 
sensitive to many cosmological parameters such as the density of barionic 
and dark matter, the Hubble constant and the index of primordial spectral 
fluctuations. Moreover, the effect of changing this parameters is similar to 
a change in a. Even though the bounds obtained using this method are not 
very stringent, it is important to consider a bound of a in the early universe. 

Martins et al }47| I48j have performed an estimation of the fine structure 
constant toghether with other cosmological parameters with the first year 
data of the Wilkinson Anisotropy Probe (WMAP) [Hlj. They established 
the following bound: 

— = 0.025 ± 0.035 (14) 
a 

3 Theoretical Model 

In this section, we review the theoretical model proposed by Damour and 
Polyakov. We derive the expression for the observable quantity — , includ- 
ing the renormalization group correction which was not considered in the 
original paper. 
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3.1 The Gravity-dilaton-matter Model 

The existence of a massless scalar field coupled to gravity and matter presents 
a host of effects that contradicts several experiments. Universal couplings 
violate the strong equivalence principle while non universal ones violate the 
universality of free fall for different composite probes. The original proposal 
of ref [U was a mechanism by which the dilaton, while still massless, is 
atracted towards some maximum of its coupling functions to the other fields 
(denoted as B a ). It is known that string loop effects associated with world- 
sheets of arbitrary genus in intermediate string sates can make B a = B a (§) 
to depend on $ in a non-monotonic way. If the extrema of this functions 
differ from each other, then the dilaton can still at present be sensitively 
changing in time, producing unacceptable time dependences of fundamental 
constants and violation to the weak equivalence principle. Conversely, if 
the coupling is universal, then the extremum will be universal as well and 
the cosmological evolution at all stages will conspire to drive the dilaton 
almost precisely towards such a value, leaving a very tiny remanent time 
dependence and very small violations of universality of free fall. 

Let us briefly describe the model. In the string frame metric denoted as 
g^v the dilaton couples in a universal multiplicative way to all other fields at 
the string tree level. However, when taking the full string-loop expansion, 
the effective action takes the general form: 



where R is the Ricci scalar corresponding to the string metric, <3? the dilaton, 



F 2 = F^F a ^ u with F^ v = d^A^-d^A^+f^A^Al includes all the gauge 



fields, and D is the gauge covariant derivative. 

Following ref.[l] we introduce the Einstein metric g^ v = CBg{^)g^ v and 
a convenient ^-dependent rescaling: 



S 




(15) 




(16) 



where B' g = dB g /d$ and C is a constant such as CB g (&o) = 1 today. Also 
we rescale the Dirac fields: 
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finally obtaining the action, decomposed into a gravitational sector (g^v,(fi) 
and a matter sector (if), A, . . .): 

A, ...} =S g ,* + S m (18) 

S m = y d 4 x^ (-i>M -^B F ($)F 2 - i$m f $) (20) 

where q = 4-7rG and G is Newton's constant (the action does not include de 
poorly-known Higgs sector). We see that the unified gauge coupling constant 
is given by: 

g-i = a^jT = kB F ($), (21) 

a function of the dilaton. Rescaling the metric means rescaling inversely 
all energies, so the string cutoff mass scale becomes dilaton dependent in 
Einstein units: 

A s ( ( j ) ) = C- 1 ^B g 1 / 2 ^)A s . (22) 

We still need to land on observational energy scales. In the case of an 
asymptotically free theory (e.g. QCD), at the one loop level the infrared 
confinement mass scale A con f is related to the cutoff scale as 

A conf oc A s exp(-8^ 2 6- 1 5 - 2 ) = C' 1 / 2 B^ 2 {((>) exp[-tof l b- 1 kB F (<[>)]A a , 

(23) 

with b dependent on the particular gauge and matter fields. 

It is assumed that the mass of any particle A depends in a non-trivial 
way on the VEV of the dilaton through the -Bs: 

m A {4>) = m A [B g (4>),B F {4>), ...] (24) 

the key hypothesis of ref.jl] being that the functions mA{4>) all have the 
same minimun 4> m that monotonically drives the dilaton to that value, as 
mentioned above. 

Another key assumption, inspired by Eg 1231 and the chiral limit in QCD 
(all energies proportional to Aqcd), is that any mass (composite or not) 
have the following form: 

m A {(t>) = ii A B- x / 2 {4>) exp[-87r 2 i/ A S((/))]A s , (25) 
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where it has been assumed that all couplings Bi depend on a common func- 
tion B through functions fi (sufficient condition to have a common ex- 
tremum); Vj\ and /j,a are pure numbers of order unity. From here we see 
that a minimum of m^(0) corresponds to a maximum of B ((/)), and hence a 
minimum of mB~ l ((j)). Following [I] we expand this function up to second 
order around the minimum: 

lnB- 1 ^) = InB- 1 ^) + - 4> m ) 2 . (26) 

From Eg 1211 we see that InoecuT °c InB^ 1 . In the next section we will 
study the cosmological model rendering the time dependence of (p. Now we 
focus on the relation between ocgut and a. First it is convenient to work at 
the unification scale Equt- ol can be written in terms of the U(l) (^) SU(2) 
coupling constants as: 

1 5 1 1 7 1 

+ — — r = (27) 



oi{Egut) 2q>i(Egut) o>2(Egut) 2o>gut 

where we have used that a\{GUT) = a2(GUT) = aour, the unification 
value. 

The renormalization group equation for the observable a is: 

1 " 7 1 +bln(^) (28) 



a(E ) 2acuT \ E J 

where b depends on the unification model, and Eo is the energy at which 
we measure a. Given that we expect Aa/a < 1CP 5 , we can write: 

Aa ( 1 1\ ( 1 1\ 

= a =i a • (29) 

a \ao a J \«o a J 

From EqE21 and assuming that Equt — A g , we expect that: 

^2EL = C i/2 B i/2^ exp [ 87r ^kB F (^)}. (30) 
Eo 

Consequently we have: 

A In J = -A In B„ + Sn 2 b- l kAB F (0). (31) 

\ Eo J 2 



Combining the last three equations we obtain: 

Aa,., (l „ b\ . _ n h AB 



a (t) = I - + 8tt^ I a Aa^ T + ^o-gA (32) 
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where «o = 1/137 is the present low energy value of the fine structure 
constant. Under the hypothesis of universality, we have B g (<p) = Bp((f)) = 
B((p), and the last equation simplifies to 

Aa ( 1 nb b\ /\B , 

— (*)= I 2+ 87r 6 + 2 ) a °^B~ = - 4AlnjB (33) 

We can estimate b from Eg 12*81 fa^rrrftn) — gn)- Regarding the string one- 
loop coefficient b which relates differente mass scales, it depends on group 
parameters of order one (aswell as b) and is positive. As we are seeking 
upper bounds for the parameter A(f> from the observational upper bounds 
on Aa, we will use the conservative lower bound A = 7/2 + 6/2. From Eq l26l 
we can then write: 

^t) = ±A K A{m)-<p m ) 2 }. (34) 

In order to proceed we must solve the cosmological time evolution of the 
dilaton, which we do in the next section. 



3.2 Cosmological Solution 

We need to know the time evolution of the dilaton, which will be apparent 
after solving Einstein equations together with the dilaton equation. The 
equation derived from the action S = + S m are: 

Rp, = 28^8^ + 2q(T lxu - -g^) (35) 



□0 = - qa (36) 

where 



[IV 



a = 2g -l/2 t_S?n (38) 

are the sources. 

The several gasses (labeled by ^4) that describe the material content of 
the Universe, which we supposed weakely interacting, can be written as: 

T*" = [ ds A m A [^x a )}u^ A 5^{x - x A ) (39) 

V 9\ x ) A J 
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and consequently 

cr(x) = -L^y [ dsAOiA[(l>{Xa)]mA[4>{Xa)]5^\x - xa) 

V9(x) a J 

= Y. a A[<P{xA)]TA{x) (40) 



where 



a A{(t>) = 71 ( 41 ) 



is a measure of the strenght of the dilaton coupling to the A particles, 
u A = dx A /dsA and Ta = ~Pa + 3Pa is the A particles contribution to the 
trace of the energy-momentum tensor. 

We consider the Friedmann cosmological model, in which ds 2 = —dt 2 + 
a(t) 2 dl 2 where dl 2 = (1- Kr 2 )- l dr 2 + r 2 (d6 2 + sra9 2 d<f) 2 ) with K = 0,1, — 1. 
Thus the gravitational and dilaton field equations simplify to: 

-3- =q(p + 3P)+24> 2 
a 

3H 2 + 3^ = 2qp + 4> 2 (42) 
4> + 3H<p = qa 

where the dot denotes the time derivative (d/dt) and H = a/a. We focus on 
the flat case K = 0. Combining eqs. ^S] and working with the logarithmic 
time p = In a + cte we have: 

—^-0" + (i _ A y = ° (43) 
(3 - 4> IZ ) p 

where A = P/p. 

In Ref.4, it is analyzed both the radiation and the matter dominated 
eras, with no cosmological constant. Even thought most of the data, belong 
to the A dominated era, we assume that the expressions for the evolution of 
the dilaton field will be similar. Therefore, we concentrate on the matter era, 
leaving the case with cosmological constant to another work. After decou- 
pling, there is no interaction between matter and radiation cosmologically 
important, so we have p = p m + p r , P = P r + P m , a = -a m {(j))(p rn - 3P m ) 
with P r = y r , P m ~ and a m ((f>) = 9ln ^ . 
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The noninteracting relativistic and non relativistic components satisfy re- 
spectively p r oc cl and p m oc m m (4>)a 3 . The evolution equation can be 
finally written as: 

0" + (l-A(p,0))0' = -[l-3A(p,0)]a rn (0) (44) 



(3 - <P' 2 ) 
where 

3A(p» = [1 + Cm m (0)e p ]- 1 (45) 

In the analysis of [3] it was concluded that the radiation era is very efficient in 
atracting the dilaton very near its minimum <j) m , and so we can approximate 
m m ((p) ~ const in X(p,(p). We choose the origin of p at matter-radiation 
equality p r (p = 0) = p m {p = 0), A(0) = l/3p r (0)/(p r (0) + p m {0)) = 1/6. 
Consequently C = m" 1 and X(p) = |(1 + e p ) _1 . 
From Eqs. (j25l26fl we obtain 



a A {4>) = (3 A (<t> - <j> m ) (46) 



where 



fo = «(40.75 -ln(^)) (47) 

and we have used the values as in Ref.[l]. Even though we will take for 
m,A = 1 GeV, the mass of the nucleon, the dark matter case m = 40 GeV 
would only change logarithmically the value of j3. Neglecting <p' 2 in Eq l44l 
we obtain: 

+ (i _ A (p))0' = -[1 - 3A(p)]/3 m (0 - 4> m ) (48) 
Denoting x = e p = a / a equiva i ence we have 

(1 + x)xd 2 x (t> + (^x + 2)d x (t> + ^(3 m (0 - m ) = 0. (49) 

This is a hypergeometric-type differential equation, with the condition of 
regularity <j>{x = 0) = (j) ra d, being 4> ra d the value of 4> at the end of the 
radiation era. The parameters are 



3 

a = - — iw 
4 

a 3 • 
p = - + IW 

4 
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7 = 2 

with w = |(/3 m - 3/8) 1/2 . The solution then reads 



<f> m ) F(a,(3,j; -x) 



(50) 



with 



F m (a,/3,7; -x) 



(51) 



(prad — 4> in 

As the argument is the red-shift, for high values we can use the asymptothic 
behaviour of F m . Unless k has the unnatural value 9.2 x 1CP 3 or less, the 
sign of (3 m — 3/8 will be positive, which we will assume. Then lo is real and 
we have: 

1 1/2 

COth(27Rij) 



F m (K, t) 



■K0j((jJ 2 + jg 



exp(— —p) cos 6 



(52) 



where 0(t) = up(t) + 2wln2 + ArgT 2 with T 2 = T (2iu) /V {2iu + |) and 
p = \n(a{t) / a eq ) . Defining 



F r («) 



■'rad tym 



(53) 



the atraction factor during the radiation era, whose value is estimated to be 
1.87 x 10" 4 x k" 9 / 4 (see RefjU) and 



Acf) = 4>i - <j) ri 



(54) 



with 4>i the "initial" (at some time ti) value of the dilaton, we can write the 
solution for <fi as: 

0{t)-0m = F m {K,t)F r (K)A0 (55) 



or more explicitly: 
<j)(t) -4> m = 1-87 x 10~ 4 x k~ 9/4 



coth(27ra;) 

7TUJ(UJ 2 + i) 



1/2 , s : 



>'eq 



cos9(t)A(f> 
(56) 
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3.3 Observational Adjustment 

Although it is tempting to work up to first or second order while analyzing 
Eg 1341 this cannot be done as we consider data spanning times comparable 
to the Hubble time. This makes the Taylor expansion of the function </>(t) 
meaningless. We will be forced then to make a non linear adjustment as we 
will see below. 

We replace the solution found in the previous section into Eg 1341 obtain- 
ing (our convention for A is A/(t) = f(t) — /(to)) 

— (t) = 1.75 x 10- 8 Ak- 7 / 2 {A(/)) 2 C(k)AV(k, t), (57) 
a 

where we have defined the functions 

„. . coth(27ru;) ,„„. 

C ^ = ( L A (58) 



V(K,t) = ( ^] 1 cos 2 (6>(t)). (59) 



and 



This entails a nonlinear adjustment of the parameters k _7//2 (A</>) 2 and k. 

Finally, it is important to note that not all data were fitted to eq. 1571 
Actually the data obtained from experiments with atomic clocks provide 
stringent bounds on — . Therefore, if we fit the derived quantity — , the 
fitting procedure and results will be almost dominated by these data. Thus 
we consider the following expression for the actual data on the variation of 
a jH: 



a ( . . 3 

— = —kHq cj tan-Otto) H — 



1.87 x l(T 4 Kf 9/4 F m ( K ,t ) Acj) 



(60) 



4 Results and Discussion 



We have performed a statistical analysis working on a x 2 function to com- 
pute the best-fit parameter values and uncertainties. For the 246 data de- 
scribed in section [5J we obtain Xmin = ^57. However, we have analized the 
contribution of each data to x 2 an d found that only 6 data have an enour- 
mous contribution while the contributions of the other 240 are of order 1 or 
less. Furthermore, excluding these data, we obtain x 2 m in = 240 for 240 data 
and 2 free parameters with no change in the value of the best fit parameters 
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and errors. Therefore, we excluded these data from the statistical analysis, 
using what we call the reduced data set. The contours of the likelihood 
functions in regions of 68 %, 95 % and 99 % of confidence level are shown 
in figure ^ Unfortunately, the contours do not close, even if we increase 
the values of n. This is due to the fact that for increasing k, A<f> has to be 
increased in order to fit the reduced data set (see eqs. |S7J EOJ). However 
from figure can find the following relation between the values of k and 
A0, that provide a good fit to the reduced data set. 

- 3.4 x 1(T 6 k < Acj>< 3.4 x 1(T 6 k (61) 

On the other hand, table El shows constraints on the values of A^ for 
fixed values of k. Again, the values show how increasing values of k, allow 
increasing values of A<f> in order to fit the reduced data set. 

We conclude that the string dilaton model proposed by Damour and 
Polyakov is able to fit almost all experimental and observational data on 
time variation of the fine structure constant. However, the constraints ob- 
tained on the free parameters (considering one or two free parameters) show 
that A(p = is also consistent with the data and in this case the model 
predicts no variation of gauge coupling constants or masses of elementary 
particles. Nevertheless, it should be noted that in any grand unified theory 
the variation of the fine structure constant is conected to the variation of the 
strong coupling constant. Thus, the prediction of this model should also be 
checked toghether with the bounds on the strong coupling constant [3*51 HUj . 
Furthermore, the data considered in section |2] should be reanalized in order 
to relate the observational values with variations in both the fine structure 
constant and the strong coupling constant. This requires a more detailed 
analysis of the QCD model considered and therefore we leave it for future 
work. 

Table 5: The table shows for fixed values of k, the best fit value and the 
corresponding la error of A<p 





A(f> 


v(A</>) 


0.01 


2.4 x 10" 14 


2.3 x 10 -13 


0.1 


-1.0 x 10" 1U 


8.0 x 10 -10 


1 


-4.2 x 10~ s 


3.2 x 10 -Y 


10 


3.6 x 10 _e 


2.7 x 10~ 5 
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Figure 1: Confidence contours for the free parameters of the DP dilaton 
model 
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